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Chapter 1 
INTRODUCTION 

In order to model many of today's complex structural systems, 
finite element models containing many thousands of degrees of freedom 
are commonly generated. Often the dimensions of the model are so large 
that a classical dynamic analysis of the system is computationally 
impossible, necessitating the use of an alternate method of analysis. 
Substructure coupling is one such analysis technique employed through- 
out the aerospace industry. 

Over the past several decades a multitude of substructure 
coupling techniques for undamped structural systems have been developed 
(Ref. 7), but comparatively few authors have concerned themselves with 
the coupling of damped systems. Neglecting the velocity terms is often 
an acceptable approximation for lightly damped structures, but is 
unacceptable for actively (or passively) controlled systems or for 
systems which develop Coriolis type forces. The object of this thesis 
is to present a general substructure coupling procedure applicable to 
systems possessing general linear damping. 

Previous papers pertaining to substructure coupling of damped 
systems include those by Hasselman and Kaplan, Hale, and Chung. The 
technique developed by Hasselman and Kaplan (Ref. 14) is an extrapo- 
lation of the popular Craig-Bampton method of component mode synthesis 
for undamped systems. Accordingly, the substructure Ritz vectors used 
by Hasselman and Kaplan are selected from the set of fixed-interface 
component modes. 
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Hale's approach to the damped substructure synthesis problem 
has been to find an applicable variational principle whose Euler 
equations are the coupled system equations of motion (Ref. 11). The 
Ritz vectors employed by Hale to represent each substructure are 
produced by a variant of subspace iteration. 

The coupling procedure presented by Chung has as its basis the 
Hamiltonian description of the system (Ref. 2). Hamilton's canonical 
equations are therefore identified as the equations of motion for the 
system. The substructure Ritz vectors utilized by Chung are a trun- 
cated set of free-interface component modes augmented by a set of 
generalized residual attachment nodes. 

The substructure coupling procedure to be presented will be 
valid for systems possessing general nonproportional, even nonsymmet- 
ric, damping terms. The coupled system equations of motion will be 
derived from a variational principle, and free-interface component 
modes along with a set of attachment modes will serve as the substruc- 
ture Ritz vectors. The presentation of the method begins in Chapter 2 
with the development of the substructure equation of motion. Chapter 3 
concerns itself with the substructure interface compatibility condi- 
tions. The actual coupling procedure is presented in Chapter 4, and 
the topic of component Ritz vectors is addressed in Chapter 5. Compu- 
tational considerations are the subject of Chapter 6. Chapter 7 
contains the results of several test problems, and conclusions and 
recommendations are drawn in Chapter 8. An appendix is provided for 
those readers unfamiliar with the properties of adjoint differential 
equations, adjoint eigenproblems, and variational principles. 


Chapter 2 

DEVELOPMENT OF THE SUBSTRUCTURE EQUATION OF MOTION 

For the purposes of this thesis, we will assume that a finite 
element model of the substructure is available. The equation of motion 
for each substructure, written in the standard form, is 



[Mjx + [C]x + [K)x 

= f 

(2.1) 

where 

x = displacement vector 

[M] 

= mass matrix 


x = velocity vector 

[C] 

= damping matrix 


x = acceleration vector 

[K] 

= stiffness matrix 


The equation of motion, when written in the above form, has 
several properties worthy of note at this point. First, and perhaps 
most importantly, Eq. (2.1) does not lend itself to either the standard 
or generalized eigenproblem form (A<f> = X<J> or A<j> = Another 
important property of Eq. {2.1} is that it is a non-self-adjoint 
ordinary differential equation. It may be demonstrated that the 
differential adjoint of Eq. (2.1) is 

[M] T J> - [C] T J + [Kl T y = f* (2.2) 

★ 

where y is the adjoint displacement vector, and f is the adjoint 
force vector. Inspection of the adjoint equation of motion shows that 
even if the defining matrices { [M] , [C], [K]) are symmetric, which is 
not assumed, the adjoint operator differs from the original differen- 
tial operator because of the sign change on the first derivative 
(damping) term. 
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To cast the equation of motion into a form which leads to a 
generalized eigenproblem, a state variable substitution will be made. 
The state variable substitution has the effect of changing the n 
second-order differential equations into 2n first-order equations. 

The transformation of the equation of motion and its adjoint 
will be accomplished by finding a variational principle which has as 
its conditions for stationarity (Euler equations) the original differ- 
ential equations (Ref. 26). The reason for applying a variational 
principle to the problem at hand is the ease of introducing constraint 
conditions in a natural manner. 

Using the concepts developed in the appendix, it may be seen 
that the following bilinear functional corresponds to the variational 
principle associated with Eqs. (2.1) and (2.2): 


IT 


1 


[y T ([M]x + [C]x + [K]x) - y T f - xY ] dt (2.3) 

0 


To find the Euler equations of this functional, the first variations of 
the functional are set to zero, i.e. 

*f 

6 ir dt = 0 (2.4) 

0 


The above expression can be expanded into the following form: 


[boundary terms] ' + 
0 

t. 


<5y T ([M]x + [C]x + [K]x - f) dt 


(2.5) 


6x T UM] T y - [C] T y + [K]\ - f*) dt = 0 


0 


11 

where the boundary terms are the by-products of the integration by 
parts. 

Since variations on the time boundary are disallowed, and 
variations Sy^ and 6 x"*” are arbitrary, the Euler equations for the 
above expression are seen to be 

[M]x + [C]x + [K]x - f = 0 (2.6a) 

and T ~ T ~ t ~ ~ 

[M]'y - [C] T y + [Kjy - f =0 (2.6b) 

which are simply the equation of motion and its adjoint. 

Up to this point, the equation of motion along with its adjoint 
have been obtained as the Euler equations of a certain functional, 
it j . For reasons discussed earlier, it is desirable to convert the 
equation of motion into a state vector form. An obvious choice for a 
state variable substitution is 


v = x (2.7) 

The above equation is equivalent to 

CM] (x - y) = 0 (2.8) 

if the mass matrix is invertible, a condition which will be assumed 
(Ref. 11). The choice of the mass matrix over other non-singular 
matrices will be made obvious shortly. 

To derive the substructure state vector formulation of the 
equation of motion, the functional is modified with the state 
variable substitution and the appending of the constraint equation to 
the functional with the use of a Lagrange multiplier vector. The 


12 


modified substructure functional, to be called it^, can be written as 


[y T ([M]y + [C]x + [K]x) - y T f - xV 


+ w T [M] (x - v) ] dt 


(2.9) 


Since there are four vectors of variables (w, y, v, x) to be considered 
independent in the tt^ functional , there will, of course, be four Euler 
equations. The four Euler equations are 


[M]x - [M]v = 0 (2.10a) 
[M]v + [C]x + [K]x - f = 0 (2.10b) 
-[M] T y - [M] T w = 0 (2.10c) 
-[M] T w - [C] T y + [K] T y - f = 0 (2.10d) 


These four equations can be conveniently stacked into the 
following two matrix equations: 



[0] [M] 


• 

V 

4. 

-[M] [0] 

_[M] [C] 


• 

x 

. 

T 

_ [0] [K] _ 

[0] M t 


r . 1 

w 


-[M] T [0] 

[M] T [C] T 


• 

y 

+ 

[0] [K]" 


(2.11a) 


(2.11b) 


Equation (2.11a) is recognized as the state variable form of the 
equation of motion (Ref. 3), and Eq. (2.11b) is simply the 

corresponding differential adjoint equation. This particular form of 
the equation of motion given in Eq. (2.11a) seems to have been utilized 
as early as 1946 by Frazer, Duncan, and Collar, who derived it 
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essentially from considerations of Hamilton's canonical equations of 


motion (Ref. 8). 

An important feature of Eqs. (2.11a) and (2.11b) is the symme- 
try of the substructure state matrices if the [M] , [C], and [K] are 
symmetric. Hale's formulation of the substructure state variable 

equation of motion (Ref. 11), although derived by a procedure similar 
to the one just presented, results in the formation of an 
unconditionally nonsymmetric substructure state matrix, which is shown 
in the appendix to be disadvantageous. 

If the equation of motion is contrasted with its adjoint, i.e. 
Eqs. (2.11a) and (2.11b), it is seen that the vector of Lagrange 
multipliers w plays the role of the adjoint state velocity. This 
idea will be utilized when compatibility between substructures is 
considered. 

Writing Eqs. (2.11a) and (2.11b) in a more compact form, we 

have 

AX + BX = F (2.12a) 

-A T Y + B T Y = F* (2.12b) 

where 


A = 

1 1 

o 

LJ 

[M] ' 

B = 

-CM] 

[0]' 


_[M] 

CCD 


. [0] 

[K] 


(2.13) 


X = 

> * 

V 

Y = 

w 


X 


X. 
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As a conclusion to this chapter, the purpose of the mass matrix 
in Eq. (2.8) v/i 1 1 be explored. The mass matrix appears in the state 
vector equation of motion, and therefore its adjoint, in two locations 
as a direct result of Eq. (2.8). Referring to the definitions of A 
and B given in Eq. (2.13), the [M] , due to its placement in 
Eq. (2.8), appears in the upper right quadrant of A , and in the upper 
left quadrant of B . If the mass matrix used in equation (2.8) is 
replaced by some other non-singular matrix, the possibility of symmetry 
in A is destroyed - a condition to be avoided if possible. 



Chapter 3 

SYSTEM CONSTRAINTS 


3.1 


Geometric Compatibility 

Consider the adjacent substructures. 



a and 6 . 


6 




The idea of interface compatibility leads directly to 


■> 

X 1 


x i 

x 2 

„ = 

X 2 

*3 

- 

a 

X 3 


(3.1) 


Equation (3.1) simply reflects the fact that the interface 

displacements of adjacent substructures are identical. This idea will 
be strictly enforced throughout the coupling procedure. Writing Eq. 
(3.1) in a more general way, we have 

■ ip (3-2) 

where the i denotes interface degrees of freedom. 

To obtain the interface degrees of freedom from the arbitrarily 
arranged displacement vector x , the concept of a "locator" matrix 
will be employed. This idea is represented by the following equation: 

x 1 = E.x (3.3) 

~ 1 ~ 
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If there are i interface degrees of freedom and n total degrees of 
freedom on the substructure, then the dimensions of Ej will be i by 
n. Essentially each row of the locator matrix selects a particular 
interface coordinate from the substructure displacement vector; hence, 
each element along the row has a value of zero, except for the element 
in the column corresponding to the location of the interface coordinate 
in the x vector. 

If Eq. (3.3) is substituted into Eq. (3.2), the result is 

(|jx) a = (|jX) B (3.4) 

which will be referred to as the displacement compatibility equation of 
substructures a and 6 . 

It is clear that the substructure velocities are subject to 
compatibility across the interface in the same manner as the displace- 
ments. However, the necessity of enforcing the velocity constraints is 
far from clear cut, and there is at least some intuitive evidence 
supporting both sides of the issue. This evidence will now be re- 
viewed. 

The supporting case for mandatory enforcement of the velocity 
constraints is usually put forth in a heuristic manner. This position 
has as its basis the physically obvious fact that the velocities are 
compatible across the interface. Since the velocities are employed as 
coordinates in the state vector formulation, any constraint condition 
on them must be included in a viable system functional. This is 
essentially the argument utilized by Chung in Ref. 2. This argument 
seems quite logical, or at least intuitively correct, but fails to 
address several key issues. 
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When the principles of Lagrangian dynamics are applied to a 
system containing constraints on displacements (and thereby velocity 
constraints) only the displacement constraints are appended to the 
Lagrangian. The velocity constraints certainly exist, but they are not 
considered in the modified Lagrangian. This situation can be con- 
sidered somewhat analogous to the substructure coupling problem at 
hand. It would appear that the constraints on the interface velocities 
are imposed throughout the system through the state vector substitution 

v = x (3.5) 

which has already been employed in the substructure equations of 
motion. In other words, the constraint on the displacements is auto- 
matically translated into a constraint on the velocities by Eq. (3.5). 
Hughes et al. (Ref. 16) have done rather extensive analysis on finite 
element systems whose displacement and velocity fields are specified 
independently, and have concluded that the velocity fields do not have 
to be coupled from element to element. If the substructures to be 
coupled are considered "superelements," then the analogy between 
adjacent substructures and adjacent finite elements is clear. 

The above concepts form the basis of Hale's position (Ref. 11) 
that the enforcement of interface velocity compatibility is optional. 
From a pragmatic point of view, test problems show that it is not 
necessary to enforce velocity compatibility. However, this is not to 
say that there are no advantages to enforcing the velocity constraints. 
There are computational advantages, and these will be considered after 
the system equations are developed. 
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To develop the actual form of the velocity compatibility 
equations, all that is required is the time derivative of Eq. (3.4), 
i .e. 


( Il*>a = (Ii*)e (3.6) 

The state vector substitution can be applied to the above equation, 
resulting in the interface velocity compatibility equation 


<IlX)a = (liX)g 


(3.7) 


Since the state vector form of the equations of motion will be 
used almost exclusively throughout the remainder of this thesis, it is 
desirable to write the above constraint equations in an appropriate 
form. For example, Eq. (3.4) can be written as 



(3.8) 


or 


E x X 
«a ~a 


Ib x ~b 


(3.9) 


Both the displacement and the velocity constraints, i.e. Eqs. 
(3.4) and (3.7), can be combined into a single matrix equation in the 
following way: 


‘e, . o' 
*1 1 = 


V 

. 

'li ; u 2 ' 


V 

0 ! E, 

_ =5 1 

a 

X 

a 

[g i h\ 

6 

X 


(3.10) 


or 
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fa (a ■ ft *6 < 3 - U > 

An additional topic concerning displacement and velocity 
compatibility between substructures is the concept of relaxed interface 
compatibility. If the substructures are represented by a set of 
partial differential equations, then satisfying the geometric 

compatibility conditions between substructures in some "average" way is 
necessary since there are an infinite number of interface degrees of 
freedom in this representation. Meirovitch and Hale (Refs. 13, 19, 20) 
have done extensive work in this area, employing a weighted residual 
approach to satisfy compatibility. These same authors, along with 
Craig and Chang (Ref. 7), have also examined relaxing the exact 
geometric compatibility conditions for substructures modeled by 
ordinary differential equations (finite element models). One effect of 
relaxing the compatibility conditions is to degrade the accuracy of the 
computed system eigenvalues. Unfortunately, an a priori estimate of 
this degradation caused by relaxing the constraints does not exist. 
Further discussion of approximate compatibility conditions will be 
postponed until the final form of the coupled system equations of 
motion has been developed. 

3.2 Force Compatibility 

Newton's Third Law (action-reaction) provides the key when a 
relationship between the interface forces is desired. This relation- 
ship is obvious when the free-body diagrams of the two adjacent sub- 
structures are considered, i.e. 
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As shown in the free-body diagrams, the interface reaction 
forces are equal in magnitude, but opposite in direction. This is 
represented by 

4 ♦ 4 * 8 ( 3 - 12 ) 

which is the equation of substructure interface reaction force com- 
patibility. 

As a conclusion to this chapter on compatibility, it should be 
recognized that the equations of displacement, velocity, and force 
compatibility developed above are obtained primarily by physical 
interpretations of the quantities involved. Since the variational 
procedure employed results in adjoint equations of motion, compatibil- 
ity equations between these "adjoint substructures" need to be devel- 
oped. Physically, it is difficult to interpret exactly what the 
adjoint substructure equations describe, but it seems logical to assume 
the adjoint compatibility equations will correspond exactly to those of 
the physical substructures. Hence the adjoint compatibility equations 
are taken to be 
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and 


E x Y 
ta la 


la U 


(3.13) 




(3.14) 



+ 



(3.15) 


Chapter 4 

DEVELOPMENT OF A SUBSTRUCTURE COUPLING PROCEDURE 

In this chapter a general substructure coupling procedure will 
be developed. Due to the assumed presence of a velocity-dependent term 
in at least one of the substructure equations of motion, all substruc- 
tures will be represented in the state vector form (Eq. (2.11)). The 
basic coupling strategy will now be outlined. 

As in the development of the substructure equations of motion, 
a variational principle will be utilized to obtain the system equations 
of motion (Ref. 11). The interface compatibility conditions developed 
in the previous section will be appended to the functional, thereby 
insuring the satisfaction of compatibility conditions. Finally, a Ritz 
approximation to the substructure state vector, X , and its adjoint 
vector, Y , will be incorporated into the functional. The Euler equa- 
tions for this final form of the system functional will be the coupled 
system equation of motion, the appropriate constraint equations, and 
the corresponding adjoint equations. 

4.1 The System Functional 

Throughout the forthcoming development, it will be assumed that 
the system under consideration is composed of two substructures, a and 
B . The developmental details of the system functional are most easily 
demonstrated when only two substructures exist, but the extrapolation 
to systems containing many components is straightforward. 
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Before the system functional is written, the substructure 
equation of motion and its adjoint will be cast in the following form: 

AX + BX = F 1 + F (4.1) 

and 

T. T *-i * 

-A'Y + B'y = F 1 + F 4.2) 

S3 ~ ~ ~ ~ 

where the force vector has been separated into two parts, the interface 

i *i 

reaction forces, F and F , and all other external forces, F 

★ i i 

and F . It should be noted that F is simply f augmented with 

zeros to length 2n. The bilinear substructure functional 

corresponding to Eqs. (4.1) and (4.2) is 

t 

* 

*3 = (AX + BX) - Y T (F 1 + F) - X T (F* 1 + F*)] dt (4.3) 

0 

In order to keep the development general at this point, the 
geometric compatibility equation to be imposed on the functional will 
be written as 

(f X) a - (E X) B = 0 (4.4) 

Note that no superscripts appear on the E matrix in this equation. 
This is intentional since it is desirable to develop the system equa- 
tions of motion with, and without, interface velocity compatibility. 
(See Eqs. (3.9) and (3.11)). 
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The system functional is simply the sum of the individual 
substructure functionals plus the appended constraint equations. 
Hence, the system functional takes the following form: 


a $ 

"3 + "3 


(4.5) 


[?i T £<!*>a - <IS>b] + C<S!> a - W] dt 


where and o ^ are Lagrange multiplier vectors. 

ct 3 

If Eq. (4.3) is substituted for ^ and ^ in Eq. (4.5), 
the following expression for the system functional results: 


t 

[ (y t [ax + bx] - yV - v t f - xV 1 - x T E*) a 

+ (Y T [AX + BX] - yV - Y T F - xV 1 - X T F*) e (4.6) 

+ ? i T [(«) 0 - * 5/ [<!!)« - <| X ) e 3 J dt 


Upon examination of the terms in the integrand of Eq. (4.6) 
which contain the interface reaction forces, it becomes clear that only 
the interface portions of the generalized displacement and force 
vectors contribute to the inner product. Thus 


yV = (E x Y) T (EV) (4.7) 

where E x is defined in Eq. (3.9). Of course, a similar expression 
exists for the adjoint force terms. 

The purpose of isolating the interface reaction forces from the 
external forces in Eqs. (4.1) and (4.2) will now be made apparent. The 
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interface reaction force terms of the system functional can be combined 
with Eq. (4.7) to yield 


-(Y T F i ) a - (Y Vjg- (X T F* i ) a - (xV 1 ) 


-(E X Y)a ( I X ^a ‘ ^| X t)e (I^Jb - (| X X)I ( 4 - 8 ) 


- (E X X)g (E x F* i ) p 

• • 

y I « 

Since E F is equivalent to f , Eq. (4.8) can be combined with the 
interface reaction force compatibility equation (i.e. Eq. (3.12)) to 
yield 


-(Y T f\- (YV)* - (xV 1 ^- (XV 1 ), = 

-[(E X Y)^ - (| X Y) gl f 1 - C(| X X)« - (| X X)gJ f* 1 


As mentioned previously, interface displacement compatibility 
will be strictly enforced throughout the coupling procedure. This 
fact, represented by Eq. (4.4), immediately allows us to conclude that 
the right hand side of Eq. (4.9) vanishes, hence 

-(yV) - (yV). - (xV 1 ) - (xV 1 ). = 0 (4.10) 

~ ~ a ~ ~ p q, ~~ p 


Equation (4.10) can be combined with Eq. (4.6), resulting 

the following form of the system functional: 
ft 


7T 


r (Y T [AX + BX] - Y T F - X T F*)_ 

L ~ ^ ^ ^ ^ u 


o 


in 


+ 


(Y T [AX + BX] - Y T F - X T F*) R (4.11) 

^ «v» -V ~ -s# 'w ~ P 

?1 T t(EX) a - (EX) b ] * o 2 T [(EY) a - (EY),]] dt 
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Notice that the interface reaction force terms no longer exist in the 
system functional. 


4.2 Introduction of Ritz Vectors 

In order to implement a reduced-order system model, a Ritz 
vector approximation to the substructure generalized displacement field 
(i.e. state vectors X and Y) will be incorporated in the system 
functional. The form of this approximation is as follows: 


and 


X = 


N x 

E 4> . T) . 

j=l XJ 


N 


= <j> . n . 

j= i~yj yj 


<J> n 
«x ~x 


$ n 
*y ~y 


(4.12) 

(4.13) 


where $ is the matrix of Ritz vectors and n is the vector of the 
time-dependent generalized coordinates. Of course, the number of Ritz 
vectors used in the approximation must not be larger than the total 
number of substructure degrees of freedom, i.e. 


N x , N y < 2n (4.14) 

The approximate system functional is formed by substituting the 
Ritz approximations into the previous expression for the system func- 
tional. This new functional takes the form 


t 
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(rj [A$ n + B<S> n ] - rj ^ F - rj $ F ) 
~y ~y ~~x~x ~ = x~x J ~y ~y - ~x -x ~ a 


+ (r|T [A$ n + B$ n ] - ^ F - n $ F ) D 

~y «y a:»x~x a^xor -!y ^y ~ -.x ~x ~ ' B 


(4.15) 


+ ?/ ^KxSx’a- fSxHxV * a J [( “y2y>a- ( “y!y^ dt 


The actual form of the Ritz vectors will be discussed in a 
subsequent chapter. 

4.3 The System Equations 

The Euler equations of the approximate system functional are 
the system equations of motion as well as the appropriate constraint 
equations. In order to insure stationarity of it , six quantities 

must be varied - nj a , nj a , , UxB* 2l T * and ~2 T * Hence the 

approximate system functional will produce six Euler equations. These 
Euler equations are easily shown to be 


(fy * *x 2x>c. + l t&K - lX 22 

(4.16a) 

* »x ?!x>6 + 1 !x2x> 6 ■ f>6 + Xy | T >6 ?2 

(4.16b) 

? T *y 2y>c, + (fl | T *y 2y>a = EV | T >a ?1 

(4.16c) 

( =x ? T !y 2y>6 * X ! T !y Vb * X E >B + X E\ 2l 

(4. 16d) 

<| *x 2x> a - ( E !x Ox>B = 2 

(4. 16e) 

( E *y V<* ' ( | !y Vb = ~ 

(4.16f ) 
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Clearly Eqs. (4.16a), (4.16b), and (4. 16e) relate to the system 
equations of motion, whereas Eqs. (4.16c), (4.16d), and (4.16f) pertain 
to the adjoint system equations of motion. Equations (4.16a) and 
(4.16b) are equivalent to the following matrix equation: 


T . 


$ 

*ya 

0 

zt 


»a 

o 

<s> n 

*y$_ 


o 

* 

*y a 

o ' 

T 

‘la 

o 

tye 


0 

$ 

~ya 

o ' 

T 

F 

~a 

o 

$ 

-yej 


.4 


0 

A. 


k 


V a 

a u 


$ 0 
*xa « 


6 * 


~x6 


n 

~xa 


~x8 


~xa 

3x3 


+ F-E $ • E$ 0 lJ a 

L » %ya i ~»v8 J 


y8- 


(4.17) 


or in abbreviated form 
* T 


A $ n 
Sbk ~x - 

+ * 1 B, . $ n 
zy sbk ~x ~ 

(4.18) 

4 T F + 

i — i 
i 

«m 

S 

4.0 

TO 

1 1 

— 1 

1 Q 
ro 


*y - 


* ~bk ’ 

are simply the block-diagonal 

matrices 


where » 
represented in Eq. (4.17). 

At this point, the equations of motion for each substructure 
have yet to be coupled in their final form. Indeed, the system 
equation represented by Eq. (4.17) is coupled only by the unknown 
vector of Lagrange multipliers (o^), which can be directly related to 


the substructure interface reaction forces. The system equation of 
motion cannot be solved until is in someway eliminated. The 
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elimination of o ^ will be achieved by applying the constraint equa- 
tions to Eq. (4.18). 


Equations (4.16e) and (4.16f) can be cast in the following 
convenient forms: 

£<!&>. <-SS,W 

and 


= 0 (4.20) 

By their very existence, the constraint equations imply that not all of 
the substructure generalized coordinates (n’s) are independent. 
If Eqs. (4.19) and (4.20) are partitioned into user-defined dependent 
and independent coordinates, the following two equations result: 


~yct| 


<i tM 


ye 


f| 

~ xa 


~ x 8 J 


(4.19) 


and 


«5&>i 


<LW 



q 





(4.21) 


(4.22) 


where the number of dependent coordinates is equal to the number of 
constraint equations to be enforced. 



Manipulation of Eq. (4.21) yields 

~xD = 1 <|!x>I Hxl 

which leads directly to the following transformation: 
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(4.23) 


?xl 


1 

Khh 

1 

2xD 


-'Ilx'o' 1 <ltx>I. 


2xl 


(4.24) 


A similar transformation of the adjoint coordinates is 


~y I 


ii 

QyD 




-yi 


(4.25) 


The transformation matrices used in Eqs. (4.24) and (4.25) will be 
designated as C and C respectively. 

A characteristic of the above transformation matrices, whose 
importance will be demonstrated shortly, is the following matrix 
identity: 




■ C <l2x>I <Itx>D ] 


il 


<I*x>o 


= 0 


(4.26) 


Naturally, an equivalent adjoint identity exists. 
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If Eq. (4.24) is substituted into Eq. (4.18), and the result 
premultiplied by C J , the following form of the system equation of 
motion results: 


£y T *y T *bk *x 3,1 + £y T !y T ? bk * x C x n xI 


Cy T t/ l * j/ [<»,), (|{,) D ] °2 


(4.27) 


In the formulation of the above equation, it has been implicitly 
assumed that the substructure Ritz vectors have been arranged in such a 
manner as to be compatible with the arrangements in Eqs. (4.21) and 
(4.22). 

The last term in Eq. (4.27) is seen to be the transpose of the 
adjoint version of Eq. (4.26); therefore, the coefficient matrix of the 
unknown Lagrange multiplier vector vanishes. The system equation of 
motion is now allowed to take its final form, 

~sys ~ sys + =sys ~ sys ~ ~sys (4.28) 


where 


-sys ~y ?y ~bk ~x ~x 


B = C T 4> T B, . $ C 
=*sys =y »y «bk «x *x 





F 


and 
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A system adjoint equation of motion can be formed in a manner 
paralleling the development of Eq. (4.28), but in practice the adjoint 
system equation appears to be of little interest in the analysis of the 
coupled system. This is not to imply that the adjoint equations which 
appear throughout the development of the coupling procedure are mere 
mathematical by-products of the variational principles employed. 
Indeed, the use of adjoint operators is central to the concept of a 
variational principle for non-self-adjoint operators such as those 
which describe the motion of generally damped systems (Ref. 9). 



Chapter 5 

COMPONENT RITZ VECTORS 


As mentioned in the previous chapter, the coupling procedure 
being developed involves a Ritz approximation to the substructure 
displacement and velocity fields. The only requirements that the Ritz 
vectors must satisfy are that they be independent and that they satisfy 
the kinematic boundary conditions of the substructure. Of course, the 
more accurately the Ritz vectors approximate the actual motion of the 
substructure, the more accurate the final coupled system will be 
modeled. 

Throughout the evolution of substructure coupling techniques, 
the formation and selection of component Ritz vectors has been a topic 
of much investigation, and justifiably so. The less computational 
effort spent defining the component Ritz vectors, the "cheaper" the 
overall coupling procedure. The vast majority of substructure coupling 
techniques employ a truncated set of component modes along with a set 
of static displacement vectors as the component Ritz vectors, hence the 
term "component mode synthesis." This is not to say that opponents of 
the use of component modes as Ritz vectors do not exist, for they do. 
This school of thought is probably best represented by Hale and 
Meirovitch, who advocate the use of "admissible vectors" as the compo- 
nent Ritz vectors (Ref. 19). These "admissible vectors" have been 
obtained in several ways, ranging in sophistication from finite element 
shape functions to variants of subspace iteration. Unfortunately, a 
formal proof that the use of admissible vectors in lieu of component 
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modes leads to comparable accuracy with less computational effort has 
yet to be demonstrated. 

For the coupling procedure being developed, a truncated set of 
substructure modes augmented by a set of attachment modes will be 
employed as the substructure Ritz vectors. The reason behind this 
particular choice of Ritz vectors is basically two-fold - first, the 
relative ease of obtaining these particular types of Ritz vectors, and 
second, the proven accuracy of methods employing normal and attachment 
modes for undamped systems (Ref. 7). 

It is clear from Eq. (4.28) that the standard and the adjoint 
Ritz vectors play equally important roles in the system equation of 
motion. Considerable computational effort will be saved if the two 
classes of Ritz vectors are chosen to be identical, i.e. 

= 9 y (5.1) 

Additionally, a practical consideration favoring the use of the 
standard eigenvectors as both types of Ritz vectors is the fact that 
the standard eigenvectors can be measured experimentally. The adjoint 
eigenvectors, however, are much more elusive, presently evading even a 
heuristic physical interpretation. 

The assumption implied by Eq. (5.1) is quite acceptable since 
the conditions that the standard Ritz vectors must satisfy are 
identical to the conditions required of the adjoint Ritz vectors 
(Ref. 11). Definitions of the particular Ritz vectors to be employed 
will now be presented. 
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5.1 Substructure Modes 

Previously it was noted that in order to form the substructure 
eigenproblem a state vector formulation was necessary. This state 
vector equation of motion was shown to have the form 

AX + BX = F (5.2) 

where A, B, F, and X are defined in Eq. (2.13). Since the coupling 
procedure deals exclusively with state vector equations, the term 
"substructure mode" will be equated with the eigenvectors of the 
homogeneous form of Eq. (5.2). 

When calculating the substructure eigenvectors, the substruc- 
tures are to be considered completely disjoint (i.e. totally isolated 
from one another). Hence the eigenvectors that are formed are the 
free-interface modes. These modes are obtained by substituting 
X = into the homogeneous state vector equation of motion, i.e. 

AA\p v + B to =0 (5.3a) 

Equation (5.3a) is recognized as a generalized eigenproblem, and can be 
readily solved by a number of algorithms. Generally the eigenvectors 
obtained from Eq. (5.3a) are complex, and commonly come in complex 
conjugate pairs. 

Logically, an adjoint eigenproblem corresponding to the adjoint 
equation of motion can be formulated. This eigenproblem is shown in 
the appendix to have the form 
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It should be repeated that the adjoint (left-hand) eigenvectors 
are not used as adjoint Ritz vectors, since Eq. (5.1) stipulates that 
the adjoint Ritz vectors are taken to be identical to the standard Ritz 
vectors. Accordingly, the standard eigenvectors serve as both standard 
and adjoint Ritz vectors, as mentioned previously. 

If desired, the substructure standard eigenproblem can be 
obtained from Eq. (5.3a) by multiplying through by A -1 . This results 
in 


M; ¥ + A” 1 B4 i =0 (5.4) 

which is recognized as a standard eigenproblem. Of course, if A - * 
does not exist, the substructure eigenproblem must be formulated in the 
generalized form. 

Upon inspection of the A matrix. A" 1 is seen to have ,the 

form 


[0] 

[M] 

[M] 

DC] 




[M] 

[Mr 1 

[0] 


(5.5a,b) 


The expression for A"* clearly shows that the condition that governs 
the existence of A - * is the existence of [Ml - * . 

S5J 

Although .the existence of A" 1 is necessary for the develop- 
ment of the standard eigenproblem, in practice the computation of A -1 
is rather inefficient and is usually replaced by one of several numer- 
ical procedures, primarily the Cholesky decomposition of A (Ref. 1). 

As a conclusion to the topic of substructure modes, a dis- 
cussion of substructure rigid-body modes will be presented. Naturally 
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an unrestrained substructure will possess rigid-body modes, but some 
question arises as to how these rigid-body modes manifest themselves in 
a state vector formulation. From the standard formulation of the 
substructure equation of motion (Eq. (2.1)), the substructure rigid- 
body modes can be defined by 


[K]x = 0 (5.6) 

where [K] is the singular stiffness matrix. 

If Eq. (5.6) defines the rigid-body modes with respect to the 
displacement coordinates, then the time derivative of Eq. (5.6) should 
lead to the rigid-body modes defined on the velocity coordinates 
(Ref. 2). Since x = v , 

[K]x = 0 (5.7) 

is equivalent to 

[K]v = 0 (5.8) 


It should be noted that the actual mode shapes generated by 
equations (5.6) and (5.8) are identical, hence the substructure state 
vector rigid-body mode set can be defined as 



<J) 

Srb 



o 

0 

23 

d> 

iRB _ 


(5.9) 


where <}> RB are the rigid-body mode shapes calculated from either Eq. 
(5.6) or (5.8). The structure of the ^ RB matrix adheres to the 
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concept of independently specified displacement and velocity fields 
(Ref. 16). Also, the dimensions of £ RB satisfy the intuitive notion 
that the number of rigid-body "modes" in the state vector formulation 
should be twice the number of rigid-body modes in the displacement 
formulation. 

A pertinent question at this point would be to ask why the 
rigid-body modes produced by the particular eigensolver employed to 
find the flexible modes discussed earlier should be replaced by the 
rigid-body modes defined in Eq. (5.9). In some cases, depending on 
what type of eigensolver is employed, the $ RB shapes may indeed be 
produced correctly by the eigensolver. Unfortunately, many times the 
<f> RB shapes are returned with distortion due to the very iteration 
process that created them (Ref. 21). It appears that the convergence 
tolerance which is quite acceptable for the flexible modes fails, in 
some instances, to provide accurate rigid-body mode shapes in the state 
vector form. Hence the need for an alternate formulation of the rigid- 
body mode shapes. 

5.2 Attachment Modes 

As stated previously, the coupling procedure being developed 
will employ a truncated set of substructure free- interface modes along 
with a set of attachment modes as the substructure Ritz vectors. 
Intuitively it is clear that if the entire set of substructure modes 
is identified as the set of Ritz vectors, then the entire "motion 
space" of the substructure will be spanned by the Ritz vectors. Since 
a reduced-order system model is desired, only a portion of substructure 
modes will be used; hence the motion space will not be completely 
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spanned. Restricting the motion space of the substructure has the 
effect of making the substructure appear stiffer than it really is. 
Said another way, the substructure has lost the flexibility represented 
by the discarded modes. Chung has shown by examples that if the 
substructure is represented only by a truncated set of low frequency 
modes, the accuracy of the system eigenvalues is rather poor (Ref. 2). 
It will be seen that if the low frequency modes are augmented by a set 
of static displacement vectors, to be called attachment modes, the 
accuracy of the system eigenvalues is significantly improved. This 
improvement is due to the fact that, in general, the attachment modes 
implicitly contain the high frequency modes which were truncated. When 
the attachment modes are included as Ritz vectors, part of the 
flexibility that they impart to the substructure is due to the implicit 
presence of the high frequency modes. 

In its most basic form, an attachment mode can be defined as 
the deflection shape the substructure takes on due to a unit load at an 
interface degree of freedom, i.e. 

[K] $ A = F (5.10) 

where [K] is the substructure stiffness matrix 

is the matrix whose columns are the attachment modes 
F is the matrix of unit interface forces 

If the substructure is restrained against rigid-body motion, the 
attachment modes are seen to be simply the columns of the flexibility 
matrix which correspond to the interface degrees of freedom. 
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To this point, the discussion of the attachment modes has only 
dealt with the displacement portion of the state vector. The velocity 
coordinates which make up the rest of the state vector will be set to 
zero since, as stated previously, the attachment modes are derived as 
the static responses to unit interface loads. 

The block-diagonal nature of B lends itself nicely to 
definition of the standard attachment modes defined in Eq. (5.10). By 
inspection, B _1 is seen to be 


-CM]' 1 

[0] 

[0] 

CK] 


assuming [M]" 1 and [K]“* exist. It is apparent that the column 
vectors of B" 1 which correspond to the displacement interface degrees 
of freedom have zero velocity portions. Hence, the standard state 
vector attachment modes are recognized as columns of generalized 
flexibility matrix, B~* . 

If the contribution of the kept low-frequency modes is removed 
from the attachment modes, thereby insuring independence of all Ritz’ 
vectors, the so-called "residual attachment modes" are formed. In 
other words, if the attachment modes developed from Eq. (5.10) are 
expressed in a modal expansion, only the portion due to the discarded 
high frequency modes will be utilized as residual attachment modes. 

Just as the standard attachment modes are defined as columns of 
the flexibility matrix, [K]~* , or more precisely the generalized 
flexibility matrix, B~* , the residual attachment modes can be defined 
from the generalized residual flexibility matrix. The residual 
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flexibility matrix for the non-self-adjoint problem at hand will now be 
developed. 

The generalized pseudo-static response problem for the sub- 
structure can be written as 


BX = F (5.12) 

The generalized displacement state vector can be expanded in terms of 
the standard right-hand eigenvectors of the substructure, i.e. 

2n 

* ■ Ui < 5 - 13 > 

where is simply the set of eigenvectors obtained from Eq. (5.3a). 
If Eq. (5.13) is substituted into Eq. (5.12), then the result pre- 
multiplied by the set of left-hand (adjoint) substructure eigenvectors 
yields 

Hr S = ^ l T F (5.14) 

where the adjoint eigenvectors (¥^) are obtained from Eq. (5.3b). 
The property of bi -orthonormal ity, developed in the appendix, allows 
Eq. (5.14) to be written as 




(5.15) 


or 


P 


tU n'l 


(5.16) 
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where the X's are the substructure eigenvalues. Substitution of Eq. 
(5.16) back into Eq. (5.13) produces the following form of the gener- 
alized displacement vector 


X 



(5.17) 


Upon examination of Eqs. (5.12) and (5.17), it becomes clear that 



(5.18) 


The right-hand side of Eq. (5.18) can be written as a summation, i.e. 



2n 
= Z 
i=l 





(5.19a) 


or 


D -1 "k $ 1R ^ iL 2 " 
B = z — — - + z 


i=l 


■ X i 


i=n k +1 


~i R hi 
- X i 


(5.19b) 


where n^ is the number of kept substructure modes. 

The last term in Eq. (5.19b) represents the flexibility pro- 


vided by the deleted high-frequency modes, and will be referred to as 
the generalized residual flexibility matrix. This residual flexibility 
matrix can be written in two ways, i.e. 
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JB 


-1 

res 


2n 

= £ 

is v’ 


tlR tlL 


(5.20) 


or 






i=l 





(5.21) 


The residual attachment modes can be defined in a manner 
similar to the standard attachment modes, the only difference being the 
use of § rec instead of B . Thus, 


$ 

*res 



(5.22) 


where $ is the matrix of residual attachment modes and F is the 
matrix containing the unit interface forces. 

As was done with standard attachment modes, the velocity 
portions of the $ vectors will be set to zero, a situation which 
occurs naturally if the damping is symmetric and proportional. 

Throughout the preceding discussion of attachment modes it was 
implicitly assumed that B"^ existed, but for an unrestrained sub- 

S3 

structure this will not be the case. Clearly an alternate method of 
defining attachment modes is in order for substructures containing 
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rigid-body degrees of freedom. Two types of attachment modes for unre- 
strained substructures will be discussed. 

When a unit force is applied to an interface degree of freedom 
of an unrestrained substructure, the substructure will exhibit both 
rigid-body and elastic motion. If the substructure is assumed to 
possess some sort of damping of elastic modes, the body will come to a 
pseudo-equilibrium state where the rigid-body acceleration still 
exists, but all elastic motion has ceased. Since no relative motion 
occurs in the substructure when the pseudo-equilibrium state has been 
reached, the damping terms will produce no internal forces. This 
argument is the basis for neglecting the damping term when the so- 
called "inertia-relief attachment modes" are calculated. The inertia- 
relief attachment modes were first introduced by Rubin (Ref. 22), but 
have been derived by Craig (Ref. 6) in a much cleaner manner. 

Since the inertia-relief attachment modes for undamped systems 
have been presented in numerous articles (Refs. 4, 6, 7, 22) they will 
not be developed here. Basically, the procedure used to calculate the 
inertia-rel ief shape is to subtract the D'Alembert forces from the 
original unit interface force vector, thereby creating a projection of 
the interface force vector which does not excite the rigid-body modes. 
This projected force vector is then applied to the substructure after 
it has been constrained in an appropriate statical ly-determinate 
manner. The deflection shapes that result are the inertia-relief 
attachment modes. Residual inertia-relief attachment modes can also be 
developed. 
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Another method of defining attachment modes for unrestrained 
substructures would be to simply constrain the substructure at user- 
defined degrees of freedom. The degrees of freedom are chosen such 
that rigid-body motion is prevented, allowing an equation similar to 
Eq. (5.10) to be utilized in defining the attachment modes (Ref. 6). 
If the degrees of freedom of the substructure are partitioned into 
three groups, the ' i ' -interface set, the ' r 1 -user-defined rigid-body 
set, and the 'o'-all other degrees of freedom set, Eq. (5.10) can be 
written as 
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(5.23) 


where R n - are the reactive forces applied to the ' r * set to prevent 
rigid-body motion. 

The top two rows of Eq. (5.23) yield 
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15. 24a) 


K ^ = F (5.24b) 

»a «a «a 

By the nature of its development, the K matrix used in Eq. (5.24b) 
is invertible; hence, <f> is readily obtained. The attachment modes 

a 

formed in this manner will be called restrained attachment modes. 
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The above development is similar to the one used in defining 
the inertia-relief attachment modes, except that when the projected 
force vector is applied to the system no reactive forces (R ri ) are 
necessary to prevent rigid-body motion. Of course, no matter which 
type of attachment mode is used for the unconstrained substructure, it 
must be augmented by a set of zeros in the velocity positions to become 
a generalized displacement Ritz vector of length 2n. 

As a conclusion to this discussion of attachment modes, it 
should be noted that the truncated substructure modes are not necessary 
to the formation of the attachment modes described (i.e. standard, 
residual, inertia-relief, or restrained attachment modes). Hence, only 
the substructure modes kept explicitly as Ritz vectors need to be 
obtained from the substructure eigenproblem. This is in contrast to 
the attachment modes utilized by Chung (Refs. 2, 3), which are formed 
directly from the high-frequency modes not used explicitly as Ritz 
vectors, therefore necessitating the complete solution of the sub- 
structure eigenproblem. 


Chapter 6 

COMPUTATIONAL CONSIDERATIONS 


In this chapter, several ideas introduced previously will be 
examined in greater detail. The topics to be considered include the 
order and type of the system eigenproblem, the issue of including the 
interface velocity constraints, and, finally, the computation of the 
state vector attachment modes. 

Equation (5.1) can be substituted into the coupled system 
equation of motion (Eq. (4.28)), resulting in 


where 


An + B n 
sys ~sys «sys ~sys 


-sys 


( 6 . 1 ) 



, cT * T fbk 
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and 


B = C T $ T B. , * C 
Ssys « » *bk * « 


Isys 


c T * t f 


Notice that in this form of the equation of motion, no distinction is 
made between the standard and the adjoint Ritz vectors, since they are 
taken to be identical. 

The system eigenproblem can be developed from Eq. (6.1) using 
the substitution n $ys = ^ S ys e< °* » 1 * e * 
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co 


A ^ 

»sys ~sys 


+ B ^ 
*sys ~sys 


0 


( 6 . 2 ) 


Two points, both important to the solution of the above system eigen- 
problem, will now be discussed. 

First, it should be clear that the A and B matrices 

^ y o -w j y o 

are complex since the component modes which are used in the formation 

of the system matrices are complex in nature. This fact will of course 

make it necessary that a complex eigensolver be available. 

A second very critical point to be considered when studying the 

system eigenproblem is the dimensions of the system matrices involved. 

From Eqs. (4.24) and (6.1) it is clear that the order of A and 

- 



N a + 


(6.3) 


where N a and Ng are the total number of Ritz vectors used to 
describe substructures a and 6 respectively, and N c is the number 
of interface constraints imposed upon the system. 

As mentioned in the section on constraints, N c can be equated 
with N. , 2N. , or N „ where N. is the number of physical 
interface degrees of freedom and N gv is the number of constraint 
equations utilized when the geometric compatibility conditions are 
satisfied in some average sense. Of course, N c = N.. implies that 

only interface displacement compatibility is satisfied, and when 
N c = 2N. both interface displacement and velocity compatibilities 
are enforced. If the geometric interface compatibility is enforced in 
an average sense, Eqs. (4.7) through (4.10) dictate that the interface 
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force compatibility must also be satisfied in an approximate fashion, a 

situation not considered in this thesis. 

The size of N g ^ s is certainly a factor when considering 

whether or not to enforce the velocity constraints. If N is 

sys 

desired to be as small as possible without altering N a or N & , it 
is clear that the velocity constraints should be enforced. Of course, 
if a larger system eigenproblem can be tolerated, the effort of enforc- 
ing the velocity constraints can be avoided by enforcing only the 
displacement compatibility. 

A computational point in favor of constraining the interface 
velocity coordinates is the conditioning of the A matrix that 
occurs when the velocity constraints are included in the formulation. 
This conditioning will be demonstrated on an axial rod with two sub- 
structures, i.e. 



Figure 6.1(a). One degree of freedom axial rod. 
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Figure 6.2(b). One degree of freedom substructures a and 8 . 

The state vector form of the equations of motion will be used 
to describe the substructures, even though the substructures will be 
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assumed to be undamped. The state vector equation of motion for the 
single element substructures can be written as 



The substructure equations of motion will be coupled twice, 
first with interface displacement and velocity coordinate constraints, 
and second with displacement constraints only. 

It may be shown that the coupled system equation of motion for 
the case of both displacement and velocity constraints is 



The characteristic equation for the above system equation of motion is 

w 2 +l = 0 (6.6) 


When only the displacement constraints are enforced, the 
coupled system equation of motion is 



The characteristic equation for this equation of motion is 


CO 2 + 1 


0 


(6.8) 
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Clearly the same solution is obtained whether or not the 
velocity constraints are enforced. However, there are two reasons why 
the eigenproblem developed from Eq. (6.5) might be easier to solve than 
the eigenproblem resulting from Eq. (6.7). First, Eq. (6.5) has fewer 
degrees of freedom than Eq. (6.7), and second, the A matrix in 
Eq. (6.5) is invertible, whereas A in Eq. (6.7) is singular. As 
mentioned previously, if the A matrix is noninvertible the gener- 
alized eigenproblem cannot be transformed into a standard eigenproblem, 
so in this case a standard eigensolver would be of no use in the 
solution of the system eigenproblem. 

As a conclusion to this chapter, the subject of attachment 
modes will be examined from a computational point of view. Of the four 
types of attachment modes developed in Chapter 5, all but the gener- 
alized residual attachment modes will be a set of unconditionally real 
vectors. For general damping, the generalized residual flexibility 
modes will be complex; however, for the special case of symmetric 
proportional damping, it can be shown that the generalized residual 
attachment modes are real . 

It has been previously noted that the adjoint eigenvectors are 
not used as Ritz vectors, thus apparently eliminating the need for the 
solution to the adjoint eigenproblem. Clearly, this is not the case if 
generalized residual attachment modes are employed as substructure Ritz 
vectors. As can be seen in Eqs. (5.21) and (5.22), the left-hand 
eigenvectors are necessary to the formation of the generalized flexi- 
bility matrix, and are therefore required to form the generalized 
residual attachment modes. 
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An alternative to the generalized residual attachment modes for 
restrained substructures are the standard attachment modes defined in 
Eq. (5.10). The solution of a set of linear algebraic equations is all 
that is required in the formation of this type of attachment mode. 

For unrestrained substructures, two types of attachment modes 
have been derived, namely the generalized inertia-relief attachment 
modes and the restrained attachment modes. The calculation of the 
inertia-relief attachment modes involves manipulation of the substruc- 
ture mass matrix followed by a solution of a set of algebraic equations 
(Ref. 6). As was required in the formation of the standard attachment 
modes, the calculation of the restrained attachment modes requires only 
the solution to a set of linear algebraic equations involving the 
substructure stiffness matrix (see Eq. (5.23)). 

It will be shown in the following chapter that the system 
accuracy obtained when using the rather unsophisticated standard or 
restrained attachment modes is comparable to the accuracy exhibited 
when generalized residual or inertia-relief attachment modes are used. 
This fact, coupled with the relative ease of computing standard and 
restrained attachment modes, forms the basis of a computational prefer- 
ence for the standard attachment modes when dealing with constrained 
substructures, and for restrained attachment modes when dealing with 
unconstrained substructures. 



Chapter 7 
EXAMPLE PROBLEMS 

The results of several test problems are presented in this 
chapter. For each problem the approximate system eigenvalues are 
compared with those of the exact system. 

EXAMPLE 1. 

The first example presented is a clamped-clamped beam with 18 
physical degrees of freedom. The beam is substructured as follows: 



A = E = L = p= l 
Clamped-Clamped System 


"T i — i i i 
Substructure a 



Substructure 3 


Figure 7.1 - 18 DOF Clamped-Clamped Beam 


53 


54 

The two cantilevered substructures have identical mass and stiffness 
properties, but the damping associated with each substructure is 
different. At the substructure level, the damping matrix is taken to 
be proportional to the stiffness matrix, as defined in the following 
equations: 

[C] a “ 48 
and 

= 96 

Of course, when the system is coupled, the system damping matrix is not 
proportional to the system stiffness matrix. 

Since the exact beam is represented by 18 physical degrees of 
freedom, the exact system eigenproblem will have 36 degrees of freedom 
due to the state vector representation of the equation of motion. The 


Ritz vectors which 

were 

used to 

represent the substructures are 

counted for below. 




Substructure 

a 

: 6 
_2 

pairs of complex conjugate modes 
real residual attachment modes 



14 

total 

Substructure 

6 

: 4 

_2 

pairs of complex conjugate modes 
real residual attachment modes 



10 

total 


Both displacement and velocity compatibility at the interface were 
enforced; therefore, the approximate system has a total of 20 degrees 
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of freedom. Table 7.1.1 presents a comparison of the system 
eigenvalues, where the complex eigenvalues are written as 

a) = a + i o)^ 

where -a is the modal damping coefficient and 
o) d is the damped natural frequency 

Also included for this example is a comparison case in which 
the same number of eigenmodes for each substructure were employed, but 
no attachment modes were included as Ritz vectors. The results for 
this second case are presented in Table 7.1.2. 

A comparison of Tables 7.1.1 and 7.1.2 clearly shows the rather 
impressive increase in accuracy provided by the inclusion of attachment 
modes as Ritz vectors. This parallels the results found by Chung 
(Ref. 2). When the attachment modes were included, the first 5 pairs 
of system eigenvalues exhibit excellent accuracy, with less than a 2% 
maximum error in the damping coefficient, and less than a 1% maximum 
error in the damped natural frequency. 



Eigenvalue Exact Approximate 0 

Number a w A a ai . % Error % Error 
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EXAMPLE 2. 

The second example is a free-free beam with the following 
geometry and substructuring: 


£ i I i i i ~ r - e= i i 

A=E = L=P=1 
Free-Free System 

r i - ~ - r ' -T - :t— — i 1 

h— c=.01k c= . 02k — H 

Substructure a 


r~:.. i i i . . . .1 

h-c=.02k-+-c=.01k-H 
Substructure 8 

Figure 7.2 - 22 DOF Free-Free Beam 

The mass and stiffness properties of the two free-free components are 
identical, but each has its own damping characteristics. As can be 
seen in Figure 7.2, the damping is not proportional at the substructure 
level; hence, the system is obviously not proportionally damped. 

Two types of attachment modes were applied to this problem, and 
the breakdown of the Ritz vectors for each case is given below. 
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Substructure 

a : 

6 

4 

2 

pairs of flexible modes 
rigid-body modes 
generalized inertia-relief 
attachment modes 



18 

total 

Substructure 

6 : 

4 

4 

2 

pairs of flexible modes 
rigid-body modes 
generalized inertia-relief 
attachment modes 



14 

total 

Substructure 

a ; 

6 

4 

2 

pairs of flexible modes 
rigid-body modes 
restrained substructure 
attachment modes 



18 

total 

Substructure 

8 : 

4 

4 

2 

pairs of flexible modes 
rigid-body modes 
restrained substructure 
attachment modes 



14 

total 


The results of Case 1 are presented in Table 7.2.1 and those of 
Case 2 in Table 7.2.2. Displacement and velocity compatibility were 
both enforced at the interface, resulting in an eigenproblem of order 
28 for both cases. The exact eigenproblem is of order 44. 

The primary goal of this example is to establish the validity 
of the definition of the rigid-body modes when a state vector approach 
is used. Also of interest is the comparison between the accuracies of 
the coupled system when inertia-relief attachment modes are used versus 
the use of restrained substructure attachment modes. 
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Although the eigenvalues corresponding to the system rigid-body 
modes were suppressed from Tables 7.2.1 and 7.2.2, the system eigen- 
solution did indeed contain 4 rigid-body modes, as would be expected 
from the state vector model of the structure. This fact, along with 
the accuracy displayed by the flexible modes, appears to justify the 
definition of the state vector rigid-body modes given in Eq. (5.9). 

When Tables 7.2.1 and 7.2.2 are compared, it is seen that both 
cases produce 8 complex conjugate pairs of system eigenvalues with a 
maximum error of approximately 5%. It should also be noted that the 
rather simple restrained attachment modes (7.2.2) produce slightly more 
accurate eigenvalues for this system than do the generalized inertia- 
relief attachment modes (Table 7.2.1). 



Eigenvalue Exact Approximate 

Number a coj a % Error % Error 
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TABLE 7.2.1 Comparison of Exact and Approximate System Eigenvalues for Example 2, 

Case 1. (Generalized inertia-relief attachment modes). All eigenvalues 
are one of a pair of complex conjugate system eigenvalues. 



Eigenvalue Exact Approximate d 

Number a &> . a a> . % Error % Error 
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TABLE 7.2.2 Comparison of Exact and Approximate System Eigenvalues for Example 2, 
Case 2. (Restrained substructure attachment modes). All eigenvalues 
are one of a pair of complex conjugate system eigenvalues. 
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EXAMPLE 3. 

The third example presented is a clamped-clamped beam with the 
geometry and substructuring shown below. 



E = A = L = P = 1 

Clamped-Clamped System 



Substructure a 


C 11 11 

Substructure 3 

Figure 7.3 - 22 DOF Clamped-Clamped Beam 

In the previous examples the damping matrix was non-propor- 
tional at the system level, but was always symmetric in form. The 
element damping matrices used in this problem are skew-symmetric in 
form, and are defined in the following equations: 

cc£ - 75 ra e 

[C]| - to tc] e 



64 


where 


[C] e 


5 

1 

1 

1 

-1 

5 

1 

1 

-1 

-1 

5 

1 

-1 

-1 

-1 

5 


The Ritz vectors utilized in the coupling procedure are now 


described. 


Substructure a : 

13 

complex modes 



2 

real standard 

attachment 




modes 



15 

total 


Substructure 6 : 

10 

complex modes 



2 

real standard 

attachment 




modes 



12 

total 



The effect of not enforcing velocity compatibility at the 
interface is examined in this example. Case 1 considers the situation 
when both displacement and velocity constraints are enforced, and Case 
2 concerns itself with displacement compatibility only. Comparisons 
between the exact and the approximate system eigenvalues for the two 
cases are presented in Tables 7.3.1 and 7.3.2 respectively. 

An interesting feature of this example, which poses no problem 
for the coupling procedure, is the fact that not all of the modes 
appear in complex conjugate pairs as in Examples 1 and 2. As can be 
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seen in either Table. 7.3.1 or 7.3.2, the first two eigenvalues have no 
oscillatory part (w d ). 

If the two tables are compared, it can be seen that both cases 
yield excellent accuracy in the low-frequency range. It should also be 
noted that, for this system, the case in which the velocity constraints 
are not enforced actually produces slightly more accurate eigenvalues 
than the case containing enforced velocity constraints. Hence, the 
fact that velocity constraints do not necessarily need to be enforced 
appears to be verified. 



Eigenvalue Exact Approximate d 

Number a aij a gj. % Error % Error 
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TABLE 7.3.1 Comparison of Exact and Approximate System Eigenvalues for Example 3, Case 



Eigenvalue Exact Approximate d 

Number a to . a to , % Error % Error 
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TABLE 7.3.2 Comparison of Exact and Approximate System Eigenvalues for Example 3, Case 
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EXAMPLE 4. 

In this final example to be presented, a pin-jointed planar 
truss structure will be considered. The geometry and substructuring of 
the system are as follows: 



Truss System 




Figure 7.4 - 36 DOF Pin-Jointed Truss. 

For the purpose of this example, only the upper and lower 
horizontal members are considered damped. The form of the element 
damping matrix for each substructure is given by the following 
equations: 
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05 [C] e 
075 [C] e 

■ i r 

-l i_ 

The nature of this problem allows us to examine the ability of 
the coupling procedure to deal with potential problems such as non- 
symmetric damping and rigid-body modes simultaneously. The substruc- 


ture Ritz vectors are given below. 



Substructure a : 

18 

complex flexible modes 


6 

real rigid-body modes 


4 

real restrained substructure 


— 

attachment modes 


28 

total 

Substructure B : 

14 

complex flexible modes 


4 

real standard attachment 




modes 


18 

total 


Both displacement and velocity coordinates are constrained at 
the interface; therefore, the approximate system eigenproblem is of 
order 38. The exact system eigenproblem is of order 72. A comparison 
between the exact and the approximate eigenvalues is given in Table 


[C] e - . 

a 



where 

[C] e = 


7.4. 
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As Table 7.4 indicates, the coupling procedure provides 8 pairs 
of complex conjugate system eigenvalues that have a maximum error of 
less than 4% in the damping coefficient. The damped natural 
frequencies are commonly two orders of magnitude larger than the 
damping coefficients, and as Table 7.4 shows, the damped natural 
frequencies are usually more accurate than the damping coefficients. 
The ability of the coupling procedure to approximate this system 
containing both nonsymmetric state matrices and rigid-body modes (at 
the substructure level) reinforces the belief that the method can be 
applied effectively to much larger complex structural systems. 
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Chapter 8 

CONCLUSIONS AND RECOMMENDATIONS 

Presented in this thesis is a procedure for the coupling of 
generally damped substructures for dynamic analysis. No assumptions of 
proportionality or even symmetry of the defining matrices (M, C, or K) 
are necessary. The method is developed from a variational principle 
which essentially results in a state vector representation of the 
system. Since low-frequency substructure modes along with a set of 
generalized attachment modes are employed as substructure Ritz vectors, 
the coupling procedure can be described as a "generalized component 
mode synthesis technique." 

The numerical test problems indicate that the coupling proce- 
dure produces accurate system eigenvalues in the low-frequency range. 
Typically, the damped natural frequencies (u^) are more accurate than 
the modal damping coefficients (a). It has also been shown numerically 
that it is not mandatory to enforce velocity compatibility between 
substructures, but that its enforcement can sometimes be advantageous. 
The rather unsophisticated attachment modes utilized appear to 
represent the truncated high-frequency modes quite well, thereby making 
the complete solution of the substructure eigenproblem unnecessary. 

On the basis of the example problems, it appears that the 
number of coupled system eigenvalues which exhibit extreme accuracy is 
at least as large as the minimum number of flexible modes used as Ritz 
vectors for either substructure. If this could be substantiated 
analytically, or even empirically, it would be a great aid to the 
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analyst, since an a priori estimate of system accuracy might dictate 
the number of substructure modes to be retained in the analysis. 

Since a knowledge of the system modes is often of great 
importance, the effect of substructuring on the accuracy of the system 
modes needs to be quantitatively explored. This suggests the estab- 
lishment of a suitable error norm which would quantify the error 
existing between the approximate and exact system eigenvectors as a 
topic for future consideration. 

Finally, any technique that efficiently produces substructure 
Ritz vectors which lead to acceptable system accuracy should be 
explored in detail. Hale's "subspace iteration" (Ref. 11) and Wilson's 
iterative procedure (Ref. 25) should both be examined further in order 
to determine their computational merits as producers of substructure 
Ritz vectors. 



Appendix 


Briefly presented in this appendix are some of the mathematical 
tools used throughout the development of the coupling procedure. The 
principal topics discussed will be adjoint differential equations, 
adjoint eigenproblems , and variational principles for non-self-adjoint 
systems. 

1. Adjoint Differential Equations 

The introduction of operator notation will prove very conve- 
nient in our discussion of adjoint differential equations, and is now 


demonstrated 

in the following example: 




Mx + Cx + Kx 

= F 

(A. 1 ) 




or 

L ( x ) = F 


(A. 2) 

where 

ii 

-in 

+ C(‘) + K( ) 

(A. 3) 


In the above equations, L is the differential operator corresponding to 
Eq. (A. 1) . The adjoint differential operator, l*, is related to L 
in the following way: 


<L(x),y> = c x,L*(y)> (A. 4) 

where < > denotes an inner product. 

It should be clear that the form of the adjoint operator is 
intimately associated with the choice of the inner product used in Eq. 
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(A. 4). The inner product usually adopted when working with differen- 


tial equations is 


<x,y> = 


y L(x) dt 


(A. 5) 


Using the above inner product, we can now find L* correspond- 
ing to the example L defined in Eq. (A. 3). This is done as follows: 


< L(x) , y > = 


y L ( x ) dt 


0 


(y^ Mx + y^ Cx + y^ Kx) dt 


[y T Mx - y T Mx + y T Cx ] * 


(A. 6) 


(y^ Mx - ^ Cx + y^ Kx) dt 


= [boundary terms] ^ 


+ 


r t 

x T (M T y - C T y + K T y) dt 

0 


Since variations on the time boundary are disallowed for this system, 
only the integral portion of the right-hand side of Eq. (A. 6) is 
non-vanishing. This fact implies that L* will be a "formal adjoint 

5 S 
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operator," i.e. an adjoint operator with no consideration of boundary 
terms. A comparison of Eqs. (A. 4) and (A. 6) easily shows that 

L*( ) = M T f) - C T f) + k t ( ) (A. 7) 

jw % m & 

Thus, the differential equation corresponding to L* is 

M T y - C T y + K T y « F* (A. 8) 

where F* represents the somewhat abstract "adjoint force" vector. 
Fortunately, a determination of F* is not necessary for the purposes 
of this thesis, and is included in Eq. (A. 8) only to provide symmetry 
with Eq. (A.l). 

Logically, if the differential operator and its adjoint are 
identical (L = L*) , the differential equation will be termed self- 
adjoint. When Eqs. (A.l) and (A. 8) are compared, it is not difficult 
to extrapolate to the fact that differential equations possessing 
odd-order derivatives will always be non-self-adjoint. 

The adjoint differential equation leads, naturally, to the 
adjoint eigenproblem, which will now be considered. 

2. The Adjoint Eigenproblem. 

Consider the following first-order linear differential 
equation: 

AX + BX = 0 (A. 9) 

t 

and its adjoint 

T* T 

-a'y + b'y = 


0 


(A. 10) 
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The usual eigenproblem can be formed from Eq. (A. 9) with the substi- 
tution 

X = ^ R e U (A. 11) 

and the adjoint eigenproblem is similarly formed from Eq. (A. 10) with 

Y = ^ L e" Xt (A. 12) 

It should be noted that the exponent in Eq. (A. 12) contains a minus 
sign, whereas the exponent in Eq. (A. 11) does not. 

The usual eigenproblem takes the form 

X A<J; r + EW r = 0 (A. 13) 

and the adjoint eigenproblem is seen to be 

XA T 4». + B T *. = 0 (A. 14) 

The adjoint eigenproblem is sometimes referred to as the "left-hand" 
eigenproblem. The reason for this terminology is obvious when the 
transpose of Eq. (A. 14) is considered, i.e. 

X, l V, l T B = 9 (A. 15) 

The eigenvectors in the above equation are on the left side of the 
and B matrices, whereas in Eq. (A. 13) the eigenvectors are on the 

as 

right side. Several important properties of the left and right-hand 
eigenproblems will now be developed. 
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Property 1. The left and right eigenvalues are identical. 

Proof: Writing the characteristic determinant of Eqs. (A. 13) 

and (A. 14) we have 

DET | AA + B | = 0 (A. 16) 


and 


DET I AA T + B T | 

S3 SS 

= DET | (AA + B) T | = 0 


(A. 17) 


Since the determinant of a matrix and its transpose are 
equal, the determinants in Eqs. (A. 16) and (A. 17) will 
yield identical characteristic equations, and therefore 
identical eigenvalues. 


Property 2. Bi -orthogonal ity, i.e. 


Proof: 


^ 1 T A^ d =0 


(A. 18) 

1 J 

1 f j 


'P 1 T B’L =0 

~ L . • 

1 J 


(A. 19) 


If Eq. (A. 13) is written for the ith eigenpair, and 
Eq. (A. 15) for the jth eigenpair, we have 


AiJ$Ri + |$ R1 = 0 (A. 20) 

and 
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Premultiplying Eq. (A. 20) by and postmul ti plying 


Eq. (A. 21) by 4; Ri results in 


X AjS*Rt + *LjS?Rf - 0 

(A. 22) 

Vli«*IH + *LJ**RI ' ° 

(A. 23) 


If Eqs. (A. 22) and (A. 23) are subtracted, the following 
equation results: 

• V ?Lj 5 *R1 ■ 0 

Assuming distinct eigenvalues (A. f A.) , then Eqs. 

* J 

(A. 24) and (A. 23) provide the bi -orthogonal ity 
conditions. 



~ Li 

**Rj 

= 0 

(A. 25) 

and 

*11 

**« 

= 0 

(A. 26) 

Property 3. Bi -orthonormal ity , 

i .e. 



If 

*l T i 

**R1 

= 1 , 

(A. 27) 

then 

*l T ( 

®*R1 

= -*1 

(A. 28) 


Proof: Premultiplying Eq. (A. 20) by *1^.. results in 

Mu . fl ?Ri + ?Li ! ~Ri - 0 (fl - 29) 
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If ~ iTi ^ ~Ri 1S norma li ze d to unity, then clearly 
4 T i itm must be equal to -A.. . 

A very useful property of the adjoint eigenvectors is their 
ability to uncouple a system of differential equations. For example, 
assume we have a system of n coupled differential equations such as 

A* + BX = F (A. 30) 

Now, the eigenvector expansion of X can be written as 

X = Y R n(t) (A. 31) 

Upon substitution of Eq. (A. 31) into Eq. (A. 30), the result can be pre- 
multiplied by V I to produce 

Il T ~ !r + ~l T § ?r 2(t) = ?iZ (A- 32 ) 

If the bi -orthonormal ity conditions are applied to the above equation, 
a set of n uncoupled differential equations of the following form 
result: 

n.(t) - A. n.(t) = F (A. 33) 

Before leaving the topic of adjoint eigenproblems , one additional point 
will be made. If the A and B matrices in Eqs. (A. 13) and (A. 14) 
are symmetric, then the left (adjoint) and right (standard) eigen- 
vectors are clearly identical. This is interesting in view of the fact 
that the differential equation from which the eigenproblems originated 
is non-self-adjoint. (See Eqs. (A. 9) and (A. 10)). This fact is very 
useful if residual attachment modes are to be employed, since their 
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calculation involves a knowledge of both left and right-hand eigen- 
vectors. 


3. Variational Principles 

A brief outline of variational principles and some of their 
uses will now be presented. The discussion is not intended to be an 
in-depth development of variational principles, but is rather aimed at 
giving the unfamiliar reader a reasonable background for understanding 
the body of this thesis and from which more rigorous presentations can 
be evolved (Ref. 9). 

For practical purposes, the key to finding a variational 
principle for a particular problem is the identification of the correct 
functional. The functional is a scalar quantity, usually expressed as 
an integral, whose conditions for stationarity are the governing 
equations for the problem. The stationarity conditions are commonly 
referred to as the "Euler equations." 

The field of structural mechanics provides the following 
simple, yet very powerful example of a variational principle. If the 
functional is written in terms of displacement variables and is equated 
with the total potential energy of the structural system, then the 
Euler equations corresponding to this functional are simply the equi- 
librium equations of the structure. This "variational procedure" is 
identical to the concept of minimum potential energy at equilibrium. 

A problem of great interest as far as the coupling procedure is 
concerned can be stated as follows: Given the governing differential 
equations, find the functional whose Euler equations reproduce the 
given differential equations. As one might expect, the form of the 
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functional is dependent upon whether or not the differential equation 
is self-adjoint or non-self-adjoint. 

If the differential equation is written as 


L(x) 



(A. 34) 


and L is a self-adjoint operator, then the functional whose Euler 
equation is Eq. (A. 34) is 
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( \ x T L ( x ) - x T F) dt 

■ 0 


(A. 35) 


As an example, it will be shown that for the following self- 
adjoint equation: 


Mx + Kx = F (A. 36) 

where M and K are symmetric, Eq. (A. 35) provides the correct 
functional . 

Upon substitution of Eq. (A. 36) into Eq. (A. 35), it is seen 
that the functional takes the form 

rt 

tt = (^x^ [Mx + Kx] - x T F) dt (A. 37) 

• 0 

The Euler equation is obtained by setting the first variation of 
equal to zero, i .e. 


6tt = 0 


(A. 38) 


Application of Eq. (A. 38) to Eq. (A. 37) results in 



(l<$x T [Mx + Kx] + jx J [M 6 x + tfx] - <Sx T F) dt = 0 (A. 39) 

• 0 

After the middle term is integrated by parts, the following expression 
is obtained: 


[x T MSX - x T M5x]J + 


6x T (Mx + Kx - F) dt = 0 


(A. 40) 


Since fix'*” is arbitrary, the following condition must exist for the 
above equation to be satisfied: 


Mx + Kx - F = 0 (A. 41) 

which is, of course, the equation we started with. 

A very useful property of the variational form of the differen- 
tial equation is the way that approximations are handled in a natural 
manner. For example, if Eq. (A. 36) is to be solved subject to the 
following Ritz vector approximation: 

x = $ n (t) (A. 42) 

then this approximation should be substituted into the variational form 
of the differential equation (Eq. (A. 39)). This substitution leads to 

•t 

( 760 "^" [M$n + K$n] + n'V [M$ 6 n + K$ 6 n] 

0 (A. 43) 

- dnW F) dt = 0 


Upon integration by parts, the above equation is reduced to 
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[n^^MMn - n^^M$6n]g + 

. (A. 44) 

ft 

6n T ( - 4^F) dt = 0 

% as as ~~ ay as ss~* a: ~ 

J 0 


Clearly, for the above expression to be satisfied for all the 

following must be true: 

4 T M$n + = 9 (A. 45) 


The above equation is simply the form of the original differential 
equation subjected to the Ritz approximation. 

An important point to consider is the case where the differen- 
tial equation represented by Eq. {A. 34) is non-self-adjoint. In this 
case the functional takes the form (Ref. 9) 


TT 


’t 

0 


(y T L(x) - F T y - F* T x) dt 


(A. 46) 


★ 

where y is the adjoint displacement vector, and F is the adjoint 
force vector. As will be shown by example, the functional in Eq. 
(A. 46) produces two matrix Euler equations. One is the original 
differential equation, and the second is the adjoint differential 
equation. 

These properties will be demonstrated in the following example: 


Ax + Bx = F 


(A. 47) 


The functional for this problem is seen to be 



t 
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it = (y T [Ax + Bx] - F T y - F* T x) dt 
■ 0 

Setting the variation of n equal to zero leads to 
•t 

tt = (<5y T [Ax + Bx] + y T [ASx + B6x] - F T 6y - F* T 6x) dt = 0 
. 0 

Integrating the second term by parts yields 

t 

[y T ^ + (6y T [Ax + Bx - F] 

0 

+ <Sx T [-A T y + B T y - F*]) dt = 0 

Since and <5x are arbitrary, the Euler equations for the func- 

tional in Eq. (A. 48) are 

Ax + Bx - F = 0 (A. 50) 

-A T y + B T y - F* = 0 (A. 51) 

Of course, Eqs. (A. 50) and (A. 51) are recognized as the original 
differential equation and its adjoint. Approximate differential 
equations can be constructed for the non-self-adjoint system in the 
same manner as for the self-adjoint case. 


(A. 48) 


(A. 49) 



REFERENCES 


1. Bathe, K.J., and Wilson, E.L., Numerical Methods in Finite 

Element Analysis , Prentice-Hall, Inc., Englewood Cliffs, NJ, 
1976. 

2. Chung, Y.T., "Application and Experimental Determination of 

Substructure Coupling For Damped Structural Systems," Ph.D. 
Dissertation, The University of Texas, Austin, TX, 1982. 

3. Chung, Y.T. and Craig, R.R., "State Vector Formulation of 

Substructure Coupling For Damped Systems," Paper No. 83-0965, 
AIAA/ASME/ASCE/AHS 24th Structures, Structural Dynamics, and 
Materials Conference , Lake Tahoe, May, 1983. 

4. Craig, R.R., Jr., Structural Dynamics: An Introduction to 

Computer Methods , John Wiley and Sons, Inc., New York, 1981. 

5. Craig, R.R., Jr. and Bampton, M.C.C., "Coupling of 
Substructures for Dynamic Analysis," AIAA Journal , Vol . 6, No. 
7, July, 1968, pp. 1313-1319. 

6. Craig, R.R., Jr. and Chang, C.J., "On the Use of Attachment 

Modes in Substructure Coupling for Dynamic Analysis," Paper No. 
77-405, AIAA/ASME/ASCE/AHS 18th Structures, Structural Dynamics 
and Materials Conference, San Diego, March, 1977, pp. 89-99. 


86 


7 


87 

7. Craig, R.R., Jr. and Chang, C.J., "Substructure Coupling For 
Dynamic Analysis and Testing," NASA CR-2781, Feb., 1977. 

8. Duncan, W.J., Collar, A.R., and Frazer, R.A., Elementary 
Matrices , Macmillan Company, New York, 1946. 

9. Finlayson, B.A. , The Method of Weighted Residuals and 
Variational Principles , Academic Press, Inc., New York, 1972. 

10. Greenberg, M.D., Foundation of Applied Mathematics , 
Prentice-Hall, Inc., Englewood Cliffs, New Jersey, 1978. 

11. Hale, A.L., "On Substructure Synthesis and Its Iterative 
Improvement For Large Nonconservative Vibratory Systems," Paper 
No. 82-0772, AIAA/ASME/ASCE/AHS 23rd Structures, Structural 
Dynamics, and Materials Conference , New Orleans, May, 1982, 
pp. 582-593. 

12. Hale, A.L. and Bergman, L.A., "The Dynamic Synthesis of General 
Nonconservative Structures From Separately Identified 
Substructure Models," Paper No. 83-0879, AIAA/ASME/ASCE/AHS 
24th Structures, Structural Dynamics, and Materials Conference , 
Lake Tahoe, May, 1983. 

13. Hale, A.L. and Meirovitch, L., "A General Substructure 
Synthesis Method for the Dynamic Simulation of Complex 
Structures," Journal of Sound and Vibration , Vol . 69, 1980, 
pp. 309-326. 



88 

14. Hasselman, T.K. and Kaplan, A., "Dynamic Analysis of Large 

Systems by Complex Modes Synthesis," Journal of Dynamic 
Systems, Measurement, and Control, Trans. ASME , Vol . 96, 

Ser. G, Sept., 1974, pp. 327-333. 

15. Hintz, R.M., "Analytical Methods in Component Modal Synthesis," 
AIAA Journal , Vol. 13, No. 8, Aug., 1975, pp. 1007-1016. 

16. Hughes, T.J.R., Hilber, H.M., and Taylor, R.L., "A Reduction 
Scheme For Problems of Structural Dynamics," International 
Journal of Solids and Structures , Vol. 12, 1976, pp. 749-767. 

17. Hurty, W.C., "Dynamic Analysis of Structural Systems Using 
Component Modes," AIAA Journal , Vol. 3, No. 4, April, 1965, 
pp. 678-685. 

18. Meirovitch, L., Computational Methods in Structural Dynamics , 
Sijthoff-Noordhoff International Publishers, Alphen ann den 
Rijn, The Netherlands, 1980. 

19. Meirovitch, L. and Hale, A.L., "A General Dynamic Synthesis For 
Structures with Discrete Substructures," Paper No. 80-0798, 
AIAA/ASME/ASCE/AHS 21st Structures, Structural Dynamics, and 
Materials Conference , Seattle, May, 1980, pp. 790-800. 

20. Meirovitch, L. and Hale, A.L., "On the Substructure Synthesis 
Method," AIAA Journal , Vol. 19, No. 7, July, 1981, pp. 940-947. 



89 

21. MSC/NASTRAN, Handbook for Dynamic Analysis , ed. M.A. Gockel , 
MSC/NASTRAN Version 63, The MacNeal -Schwendler Corp., Los 
Angeles, June, 1983. 

22. Rubin, S., "Improved Component-Mode Representation for 
Structural Dynamic Analysis," AIAA Journal , Vol . 13, No. 8, 
Aug., 1975, pp. 995-1006. 

23. Shilov, G.E., Linear Algebra , Dover Publications, Inc., New 
York, 1977. 

24. Wilkinson, J.H., The Algebraic Eigenvalue Problem , Oxford 
University Press, London, 1965. 

25. Wilson, E.L., Yuan, M.W., Dicken, J.M., "Dynamic Analysis By 
Direct Superposition of Ritz Vectors," Earthquake Engineering 
and Structure Dynamics , Vol. 10, 1982, pp. 813-821. 

26. Zienkiewicz, O.C., The Finite Element Method , McGraw-Hill Book 
Company (UK) Limited, London, 1977. 



